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i ' Abstract 

In previous work we have developed a relativistic quark model of 
mesons which is consistent with all QCD constraints at zeroth and first 
order in the heavy-quark expansion. Here we obtain first-order model 
^ ■ predictions for the differential decay spectrum, the forward-backward 

asymmetry Afb and the D* polarization parameter a in the decay 
B — > D*lu. We compare these with the predictions of QCD sum rules 
at first order. The model suggests why the corrections allowed at first 
order are small, concurrently with substantial corrections at second 
order. 
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I. INTRODUCTION 



We have recently developed a relativistic model for 5-meson semileptonic 
decays B — > D^lv, in which the hadronic matrix elements are represented 
by quark loop graphs with damping factors at the Qq- meson vertices fl], §]. 
These have the form 

where is the momentum of the light quark and P and V denote pseudoscalar 
and vector mesons. These vertices together with standard quark propagators 
determine what we will call the "full model" , in which no reference is made to 
any expansion in inverse powers of heavy-quark masses. The only parameters 
of the full model are the heavy- and light-quark masses, in terms of which the 
constants A P y and Z P y are fully determined by requiring that the meson 
self-energy functions vanish and have unit slope at the physically-measured 
meson masses. 

We may expand all quantities in the model in inverse powers of heavy 
quark masses. Severe constraints on the form of the expansion follow from 
QCD using the heavy-quark effective theory (HQET) In the heavy-quark 
limit, four of the six form factors for B — ► D^Hv become equal to a single 
universal Isgur-Wise function, while the others vanish |4j]. The twelve a 
priori independent 1 / m c and 1 / mj corrections to the form factors are given 
by specific linear combinations of the Isgur-Wise function and four additional 
universal functions |5|, |6j . The actual shapes of these universal functions are 
model-dependent, but some have model-independent values at zero recoil. 
We have shown in [1] that all such constraints are satisfied by our model at 
first order. 

The HQET restrictions arising at second order (1/toq) are available in 
[0]. We have not yet demonstrated consistency with these results, as our 
model for mesons is rather unwieldy at this order. On the other hand we 
have developed a similar model for Aq baryons || . It is more tractable than 
the meson case and is fully consistent with the HQET constraints presented 
in up to and including second order. Our approach for mesons follows the 
same principles and we thus conjecture that it is also consistent at second 
order. 

As was shown in the main feature of the full model before expansion is 
the presence of rather large deviations from the heavy quark limit predictions 
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for the zero- recoil values of the form factors h + and h,A 1 , defined in eqs. (O) 
and (|l^) below. These quantities play a crucial role in the determination of 
the Kobayashi-Maskawa element V c b, as discussed in JTDJ. Both are equal to 



1 in the heavy-quark limit and neither receives corrections at order 1/ttiq 
(Luke's theorem p). The corrections we found were traced to the effects 
of hyperfme mass splitting. For the choice of quark masses in || we found 
values of 10 to 15 per cent for the corrections at all orders, with the cor- 
rections occurring purely at second order making up about 2/3 of these full 
corrections. This translates almost directly into a 10 to 15 per cent model 
dependence in 

We believe that it is currently of interest to obtain estimates of the cor- 
rections beyond first order in 1/niQ in any model, such as ours, which is 
consistent with the heavy-quark symmetries. Many of the previously popu- 
lar models O, [T2], O] cannot be used since they do not obey the symmetry 



constraints, as described in f\M. Of course the question remains as to how 
well our model resembles QCD. Perhaps the main question has to do with 
confinement; our representation of the effects of confinement is to simply ig- 
nore the imaginary parts arising from our free quark loop calculations. We 
emphasize again, though, that this procedure is consistent with Ward iden- 
tities and heavy-quark symmetries [jlj. 

Another approach has been to carefully study the structure of the 
second-order corrections in the HQET. Here the question boils down to the 
estimation of matrix elements of various operators. Sum rules have been 
applied to the first-order corrections, but the estimates of the second-order 
corrections are much less sophisticated. In |7|] some matrix elements are es- 
timated using the Isgur-Scora-Grinstein-Wise (ISGW) nonrelativistic quark 
model pT| , but this model does not include hyperfme splitting effects. In 
fact when hyperfme splitting effects are turned off in our model (i.e. when 
g — h — in (Q) below) we find order 1/rriQ contributions to h + and at 
zero recoil which are very similar to those of the ISGW model. 

Other matrix elements involve double insertions of the chromomagnetic 
moment operator, and these are neglected in |7J since the single insertions 
appearing at first order are observed to be small. But these double insertion 
matrix elements are precisely the ones which are large in our model ||, and 
they are large concurrently with small values for the single insertion matrix 
elements. We will argue below that the small size of the first-order cor- 
rections to semileptonic decay amplitudes is the exception rather than the 
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rule. The suggested mechanism that suppresses the first-order corrections 
does not apply to higher order corrections. Even at first order large correc- 
tions are possible; for example our model and sum rules agree that the heavy 
meson decay constants receive much larger first-order corrections than the 
semileptonic decay amplitudes 

Until more reliable estimates of the second-order corrections appear, our 
model gives some indication of the possible theoretical uncertainty in the ex- 
traction of V c b. One of the objects of this paper is to provide a set of model 
predictions to be compared to future data. In particular we will define a 
"first-order model" in which only the first order corrections are kept, and 
we consider quantities which are fairly insensitive to the higher order cor- 



rections. We compare these results with QCD sum rules [15 and find that 
the differences may be large enough to make the two first-order models ex- 
perimentally distinguishable. We will also see from present data that drastic 
modifications to the model cannot be tolerated. 

In the last part of the paper we shall develop more understanding of 
the higher order corrections found in the full model. We will stress the 
consequences of holding the physical meson masses fixed in the full model, 
in particular constraints on m c and m&. 

II. THE FIRST-ORDER MODEL 

The first-order model has three parameters, besides the quark masses. 
The heavy-quark limit tuq — > 00 is characterized uniquely by the ratio A/m q , 
where m q is the light quark mass and A is the common value of Apy in 
the heavy-quark limit. The other two parameters g and h characterize the 
approach to the heavy quark limit via 

Ap,v = A(l-(6 Py h + g)—) (2) 
V m Qj 

where Sp = 3 and Sy = — 1.0 The bulk of hyperfine mass splitting is contained 
in h, a positive value of which drives the pseudoscalar mass down and the 
vector mass up by amounts in the ratio of 3 to 1. But not all effects of the 
heavy-quark spin-symmetry breaking are described by h; some breaking is 
intrinsic to the relativistic quark loop and it is present when h = 0. The three 

1 The definition of g and h here differs from that in H. 
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parameters of the first order model are free parameters, and the consistency 
with heavy-quark symmetry is valid for any values of these parameters. 

For the purpose of computing physical predictions with the first order 
model, numerical values of its parameters may be estimated as follows. As 
in 0, ||, we choose m b = 4.8 GeV, m c = 1.44 GeV and m q = 250 MeV. The 
first-order model then yields expressions for the B, B*, D and D* masses 
in terms of A, g and h. These we adjust to fit the four masses, yielding an 
optimal mass spectrum when A = 667 MeV, g = —0.13, and h = 0.19. In 
this way the model is fixed once the quark masses are given, and we find 
that changes in the first-order results are minor for any reasonable choice of 
quark masses. 

The mass difference between the meson and the heavy quark in the heavy- 



quark limit is denoted by A [16[ ; it is directly related to A |] . With the above 
values of A and m q , we have A = 504 MeV, which coincides with the QCD 
sum rule estimate in |TJ]. In addition, the Isgur-Wise function from this 
model is numerically very similar to one given in a sum rule analysis |T 



= ( 77TT ) " ■ ( 3 ) 



UJ 



The slope at u = 1 is -1.28.0 

As mentioned above, four additional universal functions, X2, Xsi an d 
£ 3 , appear at first order || [LB]. Their values are plotted for the above param- 
eter values in Fig. 1. All are relatively small compared with the Isgur-Wise 
function, which is 1 at zero recoil and approximately 0.6 at uj = 1.5. In fact 
we can begin to see how it is that dimensionless quantities appearing at first 
order are small compared to unity. \2 and £3 are spin-symmetry violating and 
conserving, respectively, but they are both independent of the parameters g 
and h for any uj They are therefore independent of the "wave function" 
distortions described by g and h which are required to fit the physical meson 
masses. Xi is spin-symmetry conserving and therefore depends only on g 
while X3 is spin-symmetry violating and depends only on h. But both xi 
and X3 are constrained to be zero at zero recoil by Luke's Theorem, and thus 
for the physical range of uj they remain small. There is thus no quantity at 
first order which is both sensitive to "wave function" distortions and nonva- 
nishing at zero recoil. In contrast, there are such dimensionless quantities 

2 More recent sum rule analyses favor a less negative slope [H. 
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appearing at second order in the heavy-quark expansion 0. Such quantities 
can be expected to be of order unity, and they thereby give corrections which 
appear to be large when compared to the first-order corrections. 

Since \2 and £3 are independent of the parameters g and h the following 
combinations of B — > D* form factors are especially interesting: 




h Aa + {m D */m B )h A2 
tlAi 




R2 



£ ( V 1 + 

These quantities have also been stressed elsewhere [jT5] . Both Ri and R 2 are 
nearly constant across the spectrum. The main difference with QCD sum 
rules is that our £3 is much smaller. Our result R\ ~ 1.3 is still in agreement 
with QCD sum rules due to the l/m& suppression of the second term. But 
R2 ~ 1.0 differs somewhat from the value of 0.8 found in QCD sum rules 
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Indeed we find, in contrast to the Isgur-Wise function, that all four uni- 
versal functions arising at first order in our model are rather different from 
those in QCD sum rules, and that they lead to different physical predictions. 
For example, we may consider the differential B — > D* spectrum in u. In the 
heavy-quark limit this goes over to \V c i,\ 2 g{uj)^{uj) 2 where g(u) is a known 
function of the meson masses and includes short distance QCD corrections. 



(We include short-distance QCD corrections |j| in all subsequent compu- 
tations.) Dividing the spectrum by g(w) and taking the square root then 
yields |V^fe|/(a;) where f{ui) goes over to in the heavy-quark limit. The 
first-order model predicts f(u>), and we obtain |V^| by fitting to the data. 

This raises an important point. The normalization f(l) = 1 at first order 
is model-independent, so it would be possible to obtain \V c b\ is a model- 
independent way if there were dcitci cit zero recoil. But the differential spec- 
trum vanishes at zero recoil, so we require an extrapolation. This extrapo- 
lation can only be accomplished by fitting some functional form to the data; 
this functional form is model-dependent and hence so is |V^,|. This model- 
dependence will diminish only as the data improves. 
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We plot in Fig. 2 the first-order model results together with ARGUS |T8 



and CLEO (T^j data for \V c b\f(uj). We see that the shape of our curve is 
steeper than that predicted by QCD sum rules, and that the present data 
data favors the steeper curve. At first order /(l) = 1, so |V^| may simply be 
read off as the intercept of the curve. We find \ V c b\ = .042 in our model and 
\Vcb\ = .037 in QCD sum rules, a difference of 13%. This large discrepancy 
is due entirely to the first-order corrections, since we have used the same 
Isgur-Wise function for both cases. 

We next present results for two integrated quantities which directly sam- 
ple the data away from zero recoil: the D* polarization parameter a and the 
forward-backward asymmetry Afb- These are plotted versus the experimen- 
tal lepton momentum cut in Figs. 3 and 4. The present experimental values 
[jTS| , [201 , HU of A FB and a are also displayed. In both cases, the predictions of 
our first-order model differ from those of QCD sum rules; better data would 
make it possible to distinguish between them. 

III. THE FULL MODEL 



If we consider the same quark masses as above then the full model yields 
[0] h+(l) = 1.107 and fiA 1 (l) = 1.155, in contrast to their model-independent 
values of 1 at first order. The value of \V c b\ is then no longer given by the 
intercept of the full model curve shown in Fig. 2, and \ V c b\ in this case is 0.38. 
In contrast, the effects of these higher order corrections manifest themselves 
only slightly in a, as shown in Fig. 3. And as seen in Fig. 4 they cancel out 
almost completely in Afb- 

We now turn to a discussion of the sensitivity of the full model results to 
the heavy-quark masses. By fixing the physical values of the meson masses, 
the quantities Z P y and A P y of Eq. are fully determined for each meson, 
as mentioned at the outset. This points to a fundamental difference in the 
interpretation of the full model versus the heavy-quark expansion. In the 
latter, the meson masses vary as the heavy-quark mass varies. But in the 
full model the physical meson masses are held fixed. Then the dependence 
of various quantities on the heavy-quark mass has no connection with the 
standard dependence in the heavy-quark expansion. In fact in the full model 
it will only be possible to produce physical meson masses when the quark 
masses are within some range. This is a feature of any realistic model of 
QCD for fixed hadron masses. In fact, the better the model, the more the 
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quark masses should be constrained to their true values. 

We can capitalize on this fact and use the full model to constrain m c 
and nit). An easily measurable quantity which has strong dependence on 
these masses is the ratio of branching ratios B(B — > D)/B(B — > D*). We 
illustrate this by plotting this ratio in Fig. 5 as a function of m c at fixed 
m b = 4.8 GeV (with no lepton momentum cut). As the data improves this 
mass dependence will translate into a constraint on m c as a function of rrif,. 

In fact, constraints on the heavy-quark masses already arise. We consider 
both B — > D* and B — > D processes and we find that the allowed region in 
the m c -mb plane is defined by 



This allowed region is displayed in Fig. 6 for a light quark mass m q = 250 
MeV. These conditions may be understood by considering the form of the 
damping factor in ([!]). They ensure that the D* is below threshold to produce 
a free charm quark and an unphysical particle with mass A#, and that the B 
is below threshold to produce a free bottom quark and an unphysical particle 
with mass Ad- (Other conditions such as m c + A# > Mb and mj + Ao» > Mb 
are less restrictive.) As we have said, the existence of an allowed region is 
expected due to the fact that the meson masses are fixed. 

Our canonical choice m c = 1.44 GeV and m b = 4.80 GeV lies close to 
a line running down the middle of the allowed region. The corrections to 
the B — > D* form factors increase as the point (m c , m&) gets closer to the 
constraint @ involving Mo* , and the corrections to the B — > D form factors 
increase as the point gets closer to the constraint (0) involving Ad. This 
effect manifests itself as the anticorrelation in the full model predictions for 
the branching ratios, as seen in Fig. 5. As another illustration we show in 
Fig. 7 the corrections to /iAi(l) and h + (l) as functions of m c with mt, = 4.8 
GeV held fixed. We again see the anticorrelation in the corrections. 

It is of interest that for no reasonable choices of quark masses are the 
corrections to /iAi(l) and h + (l) small simultaneously. In we associated 
these corrections with hyperfine splitting effects. This may be re-phrased 
in the language of Fig. 6 with the aid of Table 1 in ||. There it is seen 
that both Md* — A# and Mb — A^ decrease when hyperfine mass splitting 
is turned off. From @ and (^) this has the effect of enlarging the allowed 
region in the m c -m\, plane, and the corrections are correspondingly reduced. 



m c + A B (m b ; m q ) > M D , 
m b + A D (m c ; m q ) > M B . 



(6) 
(7) 
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It may seem most reasonable for the corrections to /^(l) and h + (l) to be 
of the same order of magnitude, and thus both in the 10% to 15% range. A 
more precise statement awaits more experimental input, as we have described. 
The main point of this paper has been to obtain definite predictions for other 
quantities, and how these predictions fare will determine how seriously this 
model should be taken. 
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APPENDIX 



We summarize in this Appendix the definitions of all relevant quantities. 
The full four-fold decay distribution in the cascade decay B° — > D* + (— > 



vi may be written as [22 



dT(B° -> D*+(-^ Dit)t-V t ) 
dui d cos 9 d cos 9* d\ 



B{D* + -> Dtt)J2 



fi(0,6*, X ) dTi 



00) 



2tt 



du 



(8) 



9 is the polar angle of the lepton measured with respect to the D*-direction 
in the {IVg) CM system, 9* is the polar angle of the D relative to the D* in 
the D* rest frame, and x is the azimuthal angle between the two decay planes 
spanned by (Dir) and (p?t). See |22|, [L8| for diagrams. B(D* + — > Dir) is the 
branching ratio T D *+^ D7T /T D *+^ a n. The zero lepton mass approximation has 
been used. The angular functions fi are listed in (P). 



U 
L 
T 
I 
P 
A 



fi 



§(1 + cos 2 #) sin"fl* 

| sin 2 9 cos 2 9* 
-^sin 2 # sin 2 ^* cos2x 
-^sin26 l sin 29* cosx 
cos 9 sin 2 9* 

16 



H; 



| 2 + \H_ 
\H a ? 



Re(H+H*_) 
±Re(H + H* a + H„H* ) 
\H + \ 2 -\H_\ 2 
-| sin 9 sin 29* cos % \^{H+H* - H-H*) 



(9) 
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The partial helicity rates dTi/duj are given by 

^l^l 2 Ml M 2 2 v^3T(l + r 2 - 2no)Hi{u), (10) 



dTj _ G 2 F W 2 
du 487T 3 



where M 12 are the B° and £)* masses and r = M2/M1. Hi are bilinear 
expressions of the three helicity amplitudes H + , H and H a describing the 
current-induced transitions B — > £)*, and are listed in (0). A set of four form 
factors may be defined by 



D*\ c^b \B°) = jMiMzhvMe^eFvfal, (11) 



D*\ c lfllb b\B°] 

-i^M x M 2 [(u + l)h Al {uj)e* 2fl - (h A2 {uj)v ltl + h A3 (u)v 2fl ) e\ ■ Vl ] ,(12) 

where e 2 is the D* polarization vector. The three helicity amplitudes are 
given in terms of the four form factors by 



H± = v /M 1 M 2 (^ + l){/ iAl T[(^-l)/(^ + l)] 1/2 ^} 

Ho = ^!g|fc + 1} {(„ - r )h Al - (a; - 1)(^ 3 + r/^)}. (13) 

The differential decay rate for 5° — > D + + £ + 77^ is given by flTj| 



^X^' ^d-cos^)^ (14) 
dud cos 9 4 gL> 

where 

It/ 12 

= ^^^^V^lCl + r 2 - 2ru,)|ff? (a;) | 2 . (15) 
M2 is now the D mass. The standard pair of form factors are defined by 

(D\ c lfl b \W) = \ .U,.U, [h+(u)(v! + v 2 )„ + h-(u)(vi - v 2 )^} . (16) 

The amplitude is given in terms of h± by 



H ° = \/1 + t*-2tu 1(1 + r)k+ " (1 " r) ^ } ' (17) 
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The total of six form factors have the following values in the heavy quark 
limit: 

hv,Ai,A 3 ,+ — £ j h A2 _ = 0. (18) 

In all of the model results quoted we include short distance QCD corrections 
f3i{uj) |TJj according to 

= h t (ur^ cD + ^{u)i{u). (19) 

We implement the momentum cut according to the prescription set out 
in l22l]; in our notation, this reads 



— 1 < cos 9 < min(cos ^(o'jpcut), 1) (20) 

where 

at \ ±-ruJ -2p cnt /M 1 

cos 9 {(jj;p cut ) = . = = • (21) 

ryw z — 1 

When converting data into values for branching ratios, the events lost due to 
the cut must be restored. In order to facilitate this we plot in Fig. 8 the full 
model prediction for the fraction of leptons having momentum greater than 
the cut, for B — > D and B — > D*. This fraction is quite insensitive to m c . 
In uncut form, the forward-backward asymmetry in the angle 9 is given 

^-krF? < 22 » 

where 

, dYi 1 + r 2 

Yi = / dm— ; u max = — — . (23) 
Ji auj 2r 

The lepton momentum cut excludes events in the extreme backward direction 
cos# — > 1; following we remove this bias by symmetrizing the definition 
of Aps with the forward hemisphere restricted by 

— min(cos6'(u;;p cut ), 1) < cos^ < 0. (24) 
In uncut form, the D* polarization parameter is given by 

a = 2^-1. (25) 
J- u 
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FIGURE CAPTIONS 



FIG. 1: First-order model predictions for universal functions £ 3 and Xi,2,3 
for A = 667 MeV , m q = 250 MeV, g = -0.13, and h = 0.19. 

FIG. 2: Predictions for the B — > D* spectrum \V c b\f(uj) (equal to |V c b|£(cj) 
in the heavy quark limit). The respective values of \V c b\ are .042 in our first 
order model, .037 in first-order QCD sum rules and .038 in our full model. 

FIG. 3: D* polarization parameter a as a function of the experimental lower 
lepton momentum cut p cut . 

FIG. 4: Forward-backward asymmetry A FB in the angle 9, as a function of 
the lepton momentum cut. The cut is performed symmetrically as explained 
in the text. 

FIG. 5: Full model prediction for B{B — > D)/B{B — > D*) as a function of 
the charm quark mass with bottom quark mass fixed at 4.8 GeV. 

FIG. 6: Allowed region in the m c -m h plane, for light quark mass m q = 250 
MeV. 

FIG. 7: Full model results for /iAi(l) and h+(l) as functions of the charm 
quark mass with bottom quark mass fixed at 4.8 GeV. 

FIG. 8: Full model prediction for fraction of leptons having momentum 
greater than a given minimum value. 
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